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Starting with an accurate pseudopotential description of the single-particle states, and following 
by configuration-interaction treatment of correlated electrons in vertically coupled, self-assembled 
InAs/GaAs quantum dot-molecules, we show how simpler, popularly-practiced approximations, de- 
pict the basic physical characteristics including the singlet-triplet splitting, degree of entanglement 
(DOE) and correlation. The mean-field-like single-configuration approaches such as Hartree-Fock 
and local spin density, lacking correlation, incorrectly identify the ground state symmetry and give 
inaccurate values for the singlet-triplet splitting and the DOE. The Hubbard model gives qual- 
itatively correct results for the ground state symmetry and singlet-triplet splitting, but produces 
significant errors in the DOE because it ignores the fact that the strain is asymmetric even if the dots 
within a molecule are identical. Finally, the Heisenberg model gives qualitatively correct ground 
state symmetry and singlet-triplet splitting only for rather large inter-dot separations, but it greatly 
overestimates the DOE as a consequence of ignoring the electron double occupancy effect. 

PACS numbers: 03.67.Mn, 73.22.Gk, 85.35.-p 



I. INTRODUCTION 

Two verticall}ii*£ or laterally coupled quantum dots 
containing electrons, holes, or an exciton constitute the 
simplest solid structure proposed for the basic gate op- 
erations of quantum computing.^ The operating prin- 
ciple is as follows: when two dots couple to each other, 
bonding and anti-bonding "molecular orbitals" (MO) en- 
sue from the single-dot orbitals of the top (T) and 
bottom (B) dots: V(°s) = Xt{s) + Xb(s) is the cr-type 
bonding and ip(a u ) = xt(s) — Xb(s) is the cr-type an- 
tibonding state. Similarly, ipi^u) = Xt(p) + Xb(p) and 
V^s) = Xt(p) — Xb{p) are the 'V bonding and anti- 
bonding states constructed from the "p" single-dot or- 
bitals of top and bottom dots, respectively. Injection of 
two electrons into such a diatomic "dot-molecule" creates 
different spin configurations such as | cr ^ T cr ^ ) or \<j^,a^), 
depicted in Fig. ^1. In the absence of spin-orbit cou- 
pling, these two-electron states are either spin-singlet or 
spin-triplet states with energy separation Js-t- Loss 
and DiVincenzo^ proposed a "swap gate" base on a sim- 
plified model, where two localized spins have Heisenberg 
coupling, H = Js-T(t)Si ■ S2. Here Si and S2 are the 
spin-i operators for the two localized electrons. The ef- 
fective Heisenberg exchange splitting Js-r{t) is a func- 
tion of time t, which is measured as the difference in the 
energy between the spin-triplet state with the total spin 
5 = 1 and the spin-singlet state with 5 = 0. The "state 
swap time" is r ~ 1/Js—t- An accurate treatment of 
the singlet-triplet splitting Js-t and the degree of en- 
tanglement carried by the two electrons is thus of out- 
most importance for this proposed approach to quantum 
computations. 

Theoretical models, however, differ in their assessment 
of the magnitude and even the sign of the singlet-triplet 
energy difference Js-t that can be realized in a quan- 
tum dot molecule (QDM) with two electrons. Most the- 



ories have attempted to model dot molecules made of 
large (50 - 100 nm), electrostatically-confined^^ dots 
having typical single-particle electronic levels separation 
of 1 - 5 meV, with larger (or comparable) inter-electronic 
Coulomb energies J ee ~ 5 meV. The central approxima- 
tion used almost universally is that the single-particle 
physics is treated via particle-in-a-box effective-mass ap- 
proximation (EM A), where multi-band and intervally 
couplings are neglected. (In this work, we will deviate 
from this tradition, see below) Many-body treatments of 
this simplified EMA model range from phenomenological 
Hubbard® or Heisenberg^' 9 models using empirical input 
parameters, to microscopic Hartree-Fock (HF) 10 ' 11 - 12 , lo- 
cal spin densities (LSD) approximationi&I^ and configu- 
ration interaction (CI) method. 2 *^ 

The LSD-EM A 1 ^ 14 can treat easily up to a few tens of 
electrons in the quantum dot molecules, but has short- 
coming for treating strongly correlated electrons, predict- 
ing for a dot molecule loaded with two electrons that 
the triplet state is below the singlet in the weak coupling 
region^ as well as incorrectly mixing singlet (spin unpo- 
larized) and triplet (spin polarized) even in the absence 
of spin-orbit coupling. Since in mean-field approaches 
like LSD or HF, the two electrons are forced to occupy 
the same molecular orbital delocalized on both dots, the 
two-electron states are purely unentangled. 

The Restricted (R)HF method (RHF-EMA) shares 
similar failures with LSD, giving a triplet as the ground 
state at large inter-dot separation. The Unrestricted (U) 
HF— corrects some of the problems of RHF by relax- 
ing the requirement of (i) two electrons of different spins 
occupying the same spatial orbital, and (ii) the single- 
particle wavefunctions have the symmetry of the exter- 
nal confining potential. The UHF-EMA correctly give 
the singlet lower in energy than the triplet^ and can 
also predict Mott localization of the electrons in the dot- 
molecule, which breaks the many-particle symmetry. 11 
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(a) Molecular orbital configurations for the dot molecule 
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FIG. 1: Possible configurations for two electrons in two ver- 
tically coupled quantum dots, (a) Spin configurations in 
the MO basis. a 3 and a u indicate the bonding and anti- 
bonding states respectively, (b) Spin configurations in the 
dot-localized basis. "T" and "B" indicate the top and bot- 
tom dot. 



However, since in UHF, the symmetry-broken wavefunc- 
tions are only the eigenstates of the z-component of total 
spin S = si + s 2 , but not of S 2 , the UHF-EMA incor- 
rectly mixes the singlet and triplet iLll For the simple 
case of dot molecules having inversion symmetry (e.g. 
molecules made of spherical dots but not of vertical lens- 
shaped dots), assuming EM A and neglecting spin-orbit 
coupling, there is an exact symmetry. For this case, Ref. 
Il^ll^ indeed were able to project out the eigenstates of 
S 2 , yielding good spin quantum numbers and lower en- 
ergy. However, for vertically coupled lens shaped quan- 
tum dots (i.e., realistic self-assembled systems) or even 
for spherical dots, but in the presence of spin-orbit cou- 
pling, there is no exact symmetry. In this case, config- 
urations with different symmetries may couple to each 
other. To get the correct energy spectrum and many- 
body wavefunctions, a further variation has to be done 
after the projection, e.g. using the Generalized Valence 
Bond (GVB) method^ For this case and other cases a 
CI approach is needed. 

The CI-EMA has been provenSii^ to be accurate for 
treating few-electron states in large electrostatic dot 
molecules, and predicts the correct ground state. Fi- 
nally, recent Quantum Monte Carlo-EMA calculations 19 
also show that the singlet is below the triplet. 

The above discussion pertained to large (50 - 100 nm) 
electrostatic-confined dots. Recently, dot molecules have 
been fabricated22i2i from self-assembled InAs/GaAs, of- 
fering a much larger Js-t- Such dot have much smaller 
confining dimensions (height of only 2-5 nm), show- 
ing a typical spacing between electron levels of 40 - 60 
meV, smaller interelectronic Coulomb energies J ee ~ 20 
meV and exchange energies of K ee ~ 3 meV. Such sin- 
gle dots have been accurately modeled^ via atomistic 
pseudopotential theories, applied to the single-particle 
problem (including multi-band and intcrvally couplings 
as well as non-parabolicity, thus completely avoiding the 
effective mass approximation). The many-particle prob- 
lem is then described via all-bound-state configuration- 



interaction method. Here we use this methodology to 
study the singlet-triplet splitting in vertically-stacked 
self-assembled InAs/GaAs dots. We calculate first the 
singlet-triplet splitting vs inter-dot separation, finding 
the singlet to be below the triplet. We then simplify 
our model in successive steps, reducing the sophistica- 
tion with which interelectronic correlation is described 
and showing how these previously practiced approxima- 
tion a 10 ! 11 ! 12 ! 13 ! 14 lead to different values of Js-t, includ- 
ing its sign reversal. This methodology provides insight 
into the electronic processes which control the singlet- 
triplet splitting in dot-molecules. 

The remainder of the paper is arranged as follows. In 
Sec.[H]we provide technical details regarding the method- 
ology we use for the calculations. We then compare 
the singlet-triplet splitting, degree of entanglement and 
correlation of two-electron states in different levels of 
approximations in Sec. II I II Finally, we summarize in 
Sec. El 



II. METHODS 
A. Geometry and strain relaxation 

We consider a realistic dot-molecule geometry^ shown 
in Fig. [21 which has recently been used in studying ex- 
citon entanglement, 4*22, and two-electron states^ Each 
InAs dot is 12 nm wide and 2 nm tall, with one mono- 
layer InAs "wetting layer", and compressively strained 
by a GaAs matrix. Even though experimentally grown 
dot molecules often have slightly different size and com- 
position profile for each dot within the molecule, here we 
prefer to consider identical dots, so as to investigate the 
extent of symmetry-breaking due to many-body effects in 
the extreme case of identical dots. The minimum-strain 
configuration is achieved at each inter-dot separation d, 
by relaxing the positions {R njQ } of all (dot + matrix) 
atoms of type a at site n, so as to minimize the bond- 
bending and bond-stretching energy using the Valence 
Force Field (VFF) method. 2 ^ This shows that both 
dots have large and nearly constant hydrostatic strain 
inside the dots which decays rapidly outside^ However, 
even though the dots comprising the molecule are geo- 
metrically identical, the strain on the two dots is different 
since the molecule lacks inversion symmetry. In fact, we 
found that the top dot is slightly more strained than the 
bottom dot. Not surprisingly, the GaAs region between 
the two dots is more severely strained than in other parts 
of the matrix, as shown in Fig. 1 of Ref. an d as the 
two dots move apart, the strain between them decreases. 



B. Calculating the single- particle states 

The single-particle electronic energy levels and wave- 
functions are obtained by solving the Schrodinger equa- 
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FIG. 2: Geometry of the two vertically coupled quantum dot 
molecule. The inter-dot distance d is measured from wetting 
layer to wetting layer. 



tions in a pseudopotential scheme, 
1 



-V 2 + V ps (r) 



V>»(r) = ei ipi(r) , 



(1) 



where the total electron-ion potential V pa (r) is a super- 
position of local, screened atomic pseudopotentials v a (r), 



i.e. 



R„ 



The pseudopoten- 



tials used for InAs/GaAs are identical to those used in 
Ref. 113 and were tested for different systemsi 23 i 27 i 28 We 
ignored spin-orbit coupling in the InAs/GaAs quantum 
dots, since it is extremely small for electrons treated here 
(but not for holes which we do not discuss in the present 
work). Without spin-orbit coupling, the states of two 
electrons are either pure singlet, or pure triplet. How- 
ever, if a spin-orbit coupling is introduced, the singlet 
state would mix with triplet state. 

Equation JQ is solved using the "linear combination 
of Bloch bands" (LCBB) method,— where the wavefunc- 
tions ipi are expanded as, 



*w=EE c - 



,k 



n.k.. 



(2) 



n,k A 



In the above equation, {4>^\ <->(')} ar e the bulk Bloch or- 
bitals of band index n and wave vector k of material A (= 
InAs, GaAs), strained uniformly to strain e*. The depen- 
dence of the basis functions on strain makes them vari- 
ationally efficient. (Note that the potential V ps (r) itself 
also has the inhomogeneous strain dependence through 
the atomic position R„. Q .) We use for the basis set e* = 
for the (unstrained) GaAs matrix material, and an av- 
erage ? value from VFF for the strained dot material 
(InAs). For the InAs/GaAs system, we use n = 2 (in- 
cluding spin) for electron states on a 6x6x28 k-mesh. 
A single dot with the geometry of Fig|2| (base=12 nm 
and hcight=2 nm) has three bound electron states (s, pi, 



and P2) and more than 10 bound hole states. The low- 
est exciton transition in the single dot occurs at energy 
1.09 eV. For the dot molecule the resulting single-particle 
states are, in order of increasing energy, the singly de- 
generated o g and a u , (bonding and antibonding combi- 
nation of the s-like single-dot orbitals), and the doubly 
(nearly) degenerated ir u and 7r ff , originating from doubly 
(nearly) degenerate "p" orbitals (split by a few meV) in 
a single dot. Here, we use the symbols g and u to de- 
note symmetric and anti-symmetric states, even though 
in our case the single-particle wavefunction are actually 
asymmetric^ We define the difference between the re- 
spective dot molecule eigenvalues as = e(a u ) — e(c 9 ) 
and A n = e{ir g ) - e(<x u ). 



C. Calculating the many-particle states 

The Hamiltonian of interacting electrons can be writ- 
ten as, 

h = E zdii>i* + ^EE r « *l ^ , (3) 

ia ijkl a,a f 

where, ipi= a u , o~ g , n u , ir g are the single-particle energy 
levels of the i-th molecular orbital, while <r, cr'=l, 2 are 
spin indices. The ri, are the Coulomb integrals between 
molecular orbitals ipi, tpj, ipk and ipi, 



1 M 



dxdv 



e(r — r')|r — r'| 



(4) 



The ./, , = rjj and K, , 



T^j are diagonal Coulomb and 



exchange integrals respectively. The remaining terms are 
called off-diagonal or scattering terms. All Coulomb inte- 
grals are calculated numerically from atomistic wavefunc- 
tions. 30 We use a phenomenological, position-dependent 
dielectric function e(r — r') to screen the electron-electron 
interaction^ 

We solve the many-body problem of Eq.© via the 
CI method, by expanding the A^-electron wavefunc- 
tion in a set of Slater determinants, |$ ei)e2) ... !ejv ) = 
^ei^L " ' '^ejvl^o), where (j>\. creates an electron in the 
state ej . The v-t\i many-particle wavefunction is then 
the linear combinations of the determinants, 



I*, 



E 

ei,e2,-" ,ejv 



A v (e-L,e 2 , ■ ■ ■ ,ejv) |<S> ei , 



(5) 



In this paper, we only discuss the two-electron problem, 
i.e. N—2. Our calculations include all possible Slater 
determinants for the six single-particle levels. 



D. Calculating pair correlation functions and 
degree of entanglement 

We calculate in addition to the energy spectrum and 
the singlet-triplet splitting Js-t also the pair correlation 
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functions and the degrees of entanglement (DOE). The 
pair correlation function P„(r,r') for an A^-particle sys- 
tem is defined as the probability of finding an electron at 
r', given that the other electron is at r, i.e., 

Pv(r,r') = J dr 3 ---dr N \V v (r,r',r 3 ,--- : r N )\ 2 , (6) 

where, ^(ri, • • • , r/v) is the iV-particle wavefunction of 
state v. For two electrons, the pair correlation function 
is just P„(r,r') = |*„(r,r')| 2 . 

The degree of entanglement (DOE) is one of the most 
important quantities for successful quantum gate opera- 
tions. For distinguishable particles such as electron and 
hole, the DOE can be calculated from Von Neumann- 
entropy formulation^ii22i22i24 However, for indistinguish- 
able particles, there are some subtletiea 35 i 36 i 37 i 38 i 39 i 40 i 41 
for defining the DOE since it is impossible to separate the 
two identical particles. Recently, a quantum correlation 
function 3 ^ has been proposed for indistinguishable par- 
ticles using the Slater decompositions^ We adapt this 
quantum correlation function to define the DOE for in- 
distinguishable fermions as, 

S = -]Tz 2 log 2 z 2 , (7) 

u 

where, z% are Slater decomposition coefficients. The de- 
tails of deriving Eq. J7J are given is Appendix [X] We 
also show in Appendix lAl that the DOE measure Eq.Q 
reduces to the usual Von Neumann-entropy formulation 
when the two-electrons are far from each other. 



III. RESULTS 

Figure shows the bonding-antibonding splitting 
A a (d) between the molecular orbitals vs inter-dot sep- 
aration d measured from one wetting layer to the other, 
showing also the value S sp = e p — e s of the splitting be- 
tween the p and s orbital energies of a single dot (i.e. 
d — > oo). The bonding-antibonding splitting decays ap- 
proximately exponentially as A CT = 2.87 exp(— d/1.15) eV 
between d ~4 - 8 nm. The result of bonding-antibonding 
splitting includes two competing effects. On one hand, 
large interdot distance d reduces the coupling between 
the two dots; on the other hand, the strain between the 
dots is also reduced, leading to a lower tunneling bar- 
rier, thus increases coupling. The local maximum of A CT 
at d=8.5 nm is a consequence of the this competition. 
Recent experiments-2i2i show the bonding-antibonding 
splitting of about 4 meV at d=11.5 nm for vertically cou- 
pled InAs/GaAs quantum dots molecules, of similar mag- 
nitude as the value obtained here(~ 1 meV), considering 
that the measured dot molecule is larger (height/base= 4 
nm/40 nm rather than 2 nm/12 nm in our calculations) 
and possibly asymmetric. We also give in Fig. [21 the in- 
terelectronic Coulomb energy Jq of a single-dot s orbital. 
We define strong coupling region as A a w 5 sp , and weak 
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FIG. 3: The bonding-antibonding splitting A CT = e(a g ) — 
e(c«) (solid line) and singlet-triplet splitting Js-t = E( 3 T,) — 
_E( 1 Eg a - ) ) (dashed line) vs inter-dot distance d. We also show 
the single-dot s, p orbitals splitting S sp — e s — e p and the s 
orbital Coulomb interaction Jq- We define "strong coupling" 
by Aa- ~ 5 sp (< 4 nm) and "weak coupling", A CT <C S sp (> 5 
nm). 

coupling region A a -C S sp . We see in Fig. [31 strong cou- 
pling for d < 4 nm, and weak coupling for d > 5 nm. In 
the weak coupling region, the 7r levels are well above the a 
levels. We also define "strong confinement" as S sp > Jc, 
and weak confinement as the reverse inequality. Figure |2I 
shows that our dot is in the strong-confinement regime. 
In contrast, electrostatic do1i£iL& are in the weak confine- 
ment regime. 

We next discuss the two-electron states in the QDMs 
and examine several different approximations which we 
call Levels 1 - 4, by comparing the properties of the 
ground states, the singlet-triplet energy separation Js-t 
and the pair correlation functions as well as the degree 
of entanglement for each state. Starting from our most 
complete model (Level 1) and simplifying it in successive 
steps, we reduce the sophistication with which interelec- 
tronic correlation is described and show how these pre- 
viously practiced approximations lead to different values 
of Js-t (including its sign reversal) , and different degree 
of entanglement. This methodology provides insight into 
the electronic features which control singlet-triplet split- 
ting and electron-electron entanglement in dot molecules. 



A. Level-1 theory: all-bound-state configuration 
interaction 

We first study the two-electron states by solving the CI 
Eq. (J5J, using all confined molecular orbitals a g , a u and 
■Kg, 7r n , to construct the Slater determinants. This gives 
a total of 66 Slater determinants. The continuum states 
are far above the bound-state, and are thus not included 
in the CI basis. Figure 0] shows some important ma- 
trix elements, including J gg (Coulomb energy of a g MO), 
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FIG. 4: Selected Coulomb integrals for molecular orbitals. 
J gg is the self-Coulomb energy of the a g orbital and J gu is 
the Coulomb energy between a g and a u orbitals, while K gu 
is the exchange energy between a g and a u orbitals. 

Jgu (Coulomb energy between a g and a u ), and K gu (ex- 
change energy between a g and a u ). The Coulomb energy 
between a M MO, J uu is nearly identical to J gg and there- 
fore is not plotted. Diagonalizing the all-bound-state CI 
problem gives the two-particle states, shown in Fig|S^. 
We show all six E states (where both electrons occupy 
the a states) and the two lowest three-fold degenerate 
3 11 U states (where one electron occupies the a g and one 
occupies one of the 7r levels). We observe that: 

(a) The ground-state is singlet for all dot-dot 

distances. However, the character of the state is very 
different at different inter-dot separation d, which can be 
analyzed by the isospin of the stated defined as the dif- 
ference in the number of electrons occupying the bonding 
(Nb) and antibonding (Nab) states in a given CI state, 
i.e. I z = (Nb — Nab)/2, where Nb and Nab ar e ob- 
tained from Eq.©. As shown in Fig© I z (d) of the l E g a) 
state is very different at different inter-dot distances: At 
small inter-dot distance, the dominant configuration of 
the ground state is \a^„, ag (both electrons occupy bond- 
ing state and Nb=2), and I z ^l. However, in the weak 
coupling region, there is significant mixing of bonding 
o g and anti-bonding states a u , and I z is smaller than 1, 
e.g I z ~ 0.2 at d— 9.5 nm. At infinite separation, where 



e B + e„ + J gu - K gu 

eg + £„ + Jg 




\ 



the bonding and antibonding states are degenerate, one 
expects I z — ► 0. 

(b) Next to the ground state, we find in FigEK the 
three-fold degenerate triplet states 3 E U , with 5*2=1, - 
1 and 0. In the absence of spin-orbit coupling, triplet 
states will not couple to singlet states. If we include 
spin-orbit coupling, the triplet may mix with the sin- 
glet state, and the degeneracy will be lifted. At large 
inter-dot distances, the ground state singlet 1 T, g and 
triplet states 3 E„ are degenerate. The splitting of to- 
tal CI energy between ground state singlet and triplet 
J S -t = E( 3 Y>) - E^Eg) is plotted in FigEJon a loga- 
rithmic scale. As we can see, Js-t also decays approxi- 
mately exponentially between 4 and 8 nm, and can be fit- 
ted as J s ~t = 5.28exp(-d/0.965) eV. The decay length 
of 0.965 nm is shorter than the decay length 1.15 nm of 
A CT . At small inter-dot separations, Js-t ~ in Fig|3 
as expected from a simple Heitler-London models 

(c) The two excited singlet states originating from the 
occupation of a u anti-bonding states, 1 S„ and 1 E g b ^ are 
further above the 3 E U state. 

(d) The lowest 3 U U states are all triplet states. They 
are energetically very close to each other since we have 
two nearly degenerate tt u MO states. In the weak cou- 
pling region, the H u states are well above the E states, 
as a consequence of large single-particle energy difference 
e(7t u ) — e(er„). However, the n u , and 1 Eg b ' ) cross at about 
4.5 nm, where the single-particle MO level tt u is still much 
higher than a u . In this case, the Coulomb correlations 
have to be taken into account. 

In the following sections, we enquire as to possible, 
popularly practiced simplifications over the all-bound- 
states CI treatment. 



B. Level-2 theory: reduced CI in the molecular 

basis 

In Level-2 theory, we will reduce the full 66x66 CI 
problem of Level- 1 to one that includes only the a g and 
a u basis, giving a 6x6 CI problem. The six many- 
body basis states are shown in FigFJJi, l a ) = |°g > °u)> 
\b)=\alai), \c)=\al,ai), \d)=\a g ,al), |e)=|aT,^), 

I/) = I°m' °m)- I n this basis set, the CI problem is reduced 
to a 6x6 matrix eigenvalue equation, 
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(8) 



where, e g and e u are the single-particle energy levels for spin-orbit coupling, the triplet states |o) and \b) are not 
the MO's \a g ) and \a u ), respectively. In the absence of 
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coupled to any other states, as required by the total spin 
conservation, and thus they are already eigenstates. The 
rest of the matrix can be solved using the integrals calcu- 
lated from Eq.(0J. The results of the 6x6 problem were 
compared (not shown) to the all-bound-state CI results: 
We find that the E states of Level-2 theory are very close 
to those of the all-bound-state CI calculations, indicat- 
ing a small coupling between a and 7r orbitals in the 
strong confinement region. We thus do not show graph- 
ically the results of Level-2. However, since we use only 
a orbitals, the II states of Level-1 (Fig[SJi) are absent 
in Level-2 theory. Especially, the important feature of 
crossover between E and H u states at 4 and 4.5 nm is 
missing. 

C. Level-3 theory: single-configuration in the 
molecular basis 

As is well known, mean-field-like treatments such as 
RHF and LSD usually give incorrect dissociation behav- 
ior of molecules, as the correlation effects are not ade- 
quately treated. Given that RHF and LSD are widely 
used in studying QMDs ) 10 i 13 i 14 it is important to under- 
stand under which circumstance the methods will succeed 
and under which circumstance they will fail in describ- 
ing the few-electron states in a QDM. In level-3 theory, 
we thus mimic the mean-field theory by further ignoring 
the off- diagonal Coulomb integrals in Eq.® of Level-2 
theory, i.e., we assume r»£= r|™ =r™=0. This ap- 
proximation is equivalent to ignoring the coupling be- 
tween different configurations, and is thus called "single- 
configuration" (SC) approximation. At the SC level, we 
have very simple analytical solutions of the two-electron 
states, 
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The energies are plotted in FigfSjD. When comparing the 
E states of the SC approach to the all-bound-state CI 
results in Figf^, we find good agreement in the strong 
coupling region for d < 5 nm (see Fig|3J|. However, the 
SC approximation fails qualitatively at larger inter-dot 
separations in two aspects: (i) The order of singlet state 
1 Eg°^ and triplet state 3 E„ is reversed (see Fig.|3>,c). (ii) 

The 1 Eg°^ and 3 E U states fail to be degenerate at large 
interdot separation. This lack of degeneracy is also ob- 
served for 1 Eg fc ' 1 and 1 E U . These failures are due to the 
absence of correlations in the SC approximation. Indeed 
as shown in FigEl the accurate Level-1 ground state sin- 
glet has considerable mixing of anti-bonding states, i.e. 
I z — > at large d. However, in the SC approximation 
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FIG. 5: (Color online) (a) Two-electron states calculated 
from CI using all confined MO from LCBB (Level-1), includ- 
ing the singlet 1 £g , 1 E U , states and the 3-fold degen- 
erated triplet states 3 E„ as well as two 3-fold degenerated 
triplet states 3 T1 U . (b) Two electron states calculated from 
the single-configuration approximation (Level-3). (c) Com- 
parison of the singlet-triplet splitting calculated from Levels 
1, 3 and 4 theories. 



both electrons are forced to occupy the o~ g orbital in the 

lowest singlet state 1 T,g as a consequence of the lack of 
the coupling between the configuration |e) of Fig^ and 
other configurations. As a result, in Level-3 theory, the 
isospins are forced to be I z =l for 1 Eg Q ' 1 at all inter-dot 
distances d, which pushes the singlet energy higher than 
the triplet. 
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Isospin of 1 £f state 

i ' 1 1 1 1 1 1 r 



Level-3 




Inter-dot distance (nm) 



FIG. 6: Isospin, defined as the difference in the number 
of electrons occupying the bonding (Nb) and antibonding 
(Nab) states, of the 1 'Eg state in Level- 1 and Level-3 theo- 
ries. 



where, r\ — (l,p) and x\,a creates an electron in the /=(s, 
p, • • • ) orbital on the p={T, B) dot with spin a that 
has single-particle energy e v . Here, t VlTj2 is the coupling 
between the r\\ and 772 orbitals, and rj^'JK is the Coulomb 
integral of single-dot orbitals Xm , Xr n , Xvz an0 - Xm ■ 

We wish to construct a Hubbard Hamiltonian whose 
parameters are taken from the fully atomistic single- 
particle theory. To obtain such parameters in Eq.|TJ3J) 
including e v , t mV2 and T^^, we resort to a Wannier-like 
transformation, which transform the "molecular" orbitals 
(Fig^i) into single-dot "atomic" orbitals (Fig^>). The 
latter dot-centered orbitals are obtained from a unitary 
rotation of the molecular orbitals ipi, i.e, 

X v =^^,i^i , (14) 
i=i 

where, ipi is the 7-th molecular orbitals, Xn 1S t ne sm ~ 
gle dot-centered orbitals, and U are unitary matrices, i.e. 
U^U = I. We chose the unitary matrices that maximize 
the total orbital self-Coulomb energy. The procedure of 
finding these unitary matrices is described in detail in 
Appendix The dot-centered orbitals constructed this 
way are approximately invariant to the change of cou- 
pling between the dots^i Once we have the U matrices, 
we can obtain all the parameters in Eq. (|13ll by trans- 
forming them from the molecular basis. The Coulomb 
integrals in the new basis set are given by Eq. I|B2|1 . while 
other quantities including the effective single-particle lev- 
els e v for the 77-th dot-centered orbital, and the coupling 



D. Level-4 theory: Hubbard model and Heisenberg 
model in a dot-centered basis 

The Hubbard model and the Heisenberg model are of- 
ten used£ to analyze entanglement and gate operations 
for two spins qbits in a QDM. Here, we analyze the ex- 
tent to which such approaches can correctly capture the 
qualitative physics given by more sophisticated models. 
Furthermore, by doing so, we obtain the parameters of 
the models from realistic calculations. 



1. Transforming the states to a dot-centered basis 

Unlike the Level 1-3 theories, the Hubbard and the 
Heisenberg models are written in a dot-centered basis 
as shown in Fig^>, rather than in the molecular basis 
of Fig^,. In a dot-centered basis, the Hamiltonian of 
Eq.® can be rewritten as, 



(13) 

I 

between the 771-th and 772-th orbitals t mri2 can be obtained 
from, 

e v = {Xr,\T\Xn) =^2, U n,i U v,i e i > ( 15 ) 

i 

tmm = (Xm\T\Xr, 2 ) = J2 U m,i U i2,i e i . ( 16 ) 

i 

where ej is the single-particle level of the i-th molecu- 
lar orbital, and T is kinetic energy operator. Using the 
transformation of Eq.JTSJ, Eq.JTU and Eq.jB2j), we cal- 
culate all parameters of Ea. (|13fl . Figure [7^,, shows the 
effective single-dot energy of the "s" orbitals obtained 
in the Wannier representation for both top and bottom 
dots. We see that the effective single-dot energy levels 
increase rapidly for small d. Furthermore the energy lev- 
els for the top and bottom orbitals are split due to the 
strain asymmetry between the two dots. We compute the 
Coulomb energies Jtt, Jbb of the "s" orbitals on both 
top and bottom dots, and the inter-dot Coulomb and 
exchange energies Jtb and Ktb and plot these quanti- 
ties in Fig. 03. Since Jtt and Jbb are very similar, we 
plot only Jtt- As we can see, the Coulomb energies of 
the dot-centered orbitals are very close to the Coulomb 
energy of the s orbitals of a isolated single dot (dashed 
line). The inter-dot Coulomb energy Jtb has comparable 
amplitude to Jtt and decays slowly with distance, and 
remain very significant, even at large separations. How- 
ever, the exchange energy between the orbitals localized 
on top dot and bottom dot i^TB is extremely small even 
when the dots are very close. 



V2 -wt T/t y /V 
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2. "First-principles" Hubbard model and Heisenberg model: 
Level-4 

In level-4 approximation, we use only the "s" orbital in 
each dot. Figure ^Jd shows all possible many-body basis 



functions of two electrons, where top and bottom dots are 
denoted by "T" and "B" respectively. The Hamiltonian 
in this basis set is, 



H = 













et + e-B + Jte 







Ktb 






ct + e B + ./tb 

— Kjb 

+ T-1 BB 



T B 
TT 

TB 






-a'tb t - r 

et + e B + ./tb — t + r 

-t+fg^ 2e B +./B 

-t + tTI o 



TB 
Li B 
~ TB 








2eT + Jtt 



(17) 



where t = txB and to simplify the notation, we ignore the 
orbital index "s" . If we keep all the matrix elements, the 
description using the molecular basis of Fig. \Ql and the 
dot localized basis of Fig. ^ are equivalent, since they are 
connected by unitary transformations. We now simplify 
Ea. (|17|l by ignoring the off-diagonal Coulomb integrals. 
The resulting Hamiltonian is the single-band Hubbard 
model. Unlike Level-3 theory, in this case, ignoring off- 
diagonal Coulomb integrals (but keeping hopping) can 
still give qualitatively correct results, due to the fact that 
off-diagonal Coulomb integrals such as <C t, and the 
correlation is mainly carried by inter-dot hopping t. We 
can further simplify the model by assuming ex = ee = e; 
Jtt = Jbb = U; and let Jtb = V, Ktb = K. We can 
then solve the simplified eigenvalue equation analytically. 
The eigenvalues of the above Hamiltonian are (in order 
of increasing energy) : 

1. Ground state singlet 1 E 
1 



(a) Effective single-particle levels 



(a) 



E = 2e+-[U + V + K- Vl6* 2 + (U —V — K) 2 ] (18) 



2. Triplet states (three- fold degenerate) 3 S M 
E = 2e + V - K 



3. Singlet 1 S i 



4. Singlet 



E = 2e + U 



(19) 



(20) 



E = 2e + -[U + V + K+ ^16t 2 + {U -V- K) 2 ] (21) 

In the Hubbard limit where Coulomb energy U 3> £, 
the singlet-triplet splitting J S -t = E( 3 T,) - i?( 1 S g ) ~ 
4t 2 /(U — V), which reduces the model to the Heisenberg 
model 



H 



U-V 



St ■ Sb 



(22) 



where St and Sb are the spin vectors on the top and 
bottom dots. The Heisenberg model gives the correct 
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FIG. 7: (a) Effective single-particle energy levels of s orbitals 
localized on the top (er) and bottom (es) dots, (b) Intra-dot 
Coulomb energy Jtt , Jbb , inter-dot Coulomb energy Jtb and 
inter-dot exchange energy -Ktb (magnified by a factor 100). 
The dashed line gives the single-dot s orbital self- Coulomb 
energy J c . 



order for singlet and triplet states. The singlet-triplet 
splitting Js-t = 4i 2 /(£7 - V) is plotted in Fig. Et and 
compared to the results from all-bound-state CI calcula- 
tions (Level-1), and single-configuration approximations 
(Level-3). As we can see, at d > 6.5 nm, the agreement 
between the Heisenberg model with CI is good, but the 
Heisenberg model greatly overestimates Js-t at d < 6 
nm. 
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Pair correlation functions of 1 xl a) state 





Pair correlation functions of 1 £l b) state 



Level- 1 





FIG. 8: (Color online) Comparison of pair correlation func- 
tions calculated from (a) Level- 1, (b) Level-3 theory for the 
1 Hg state and (c) Level-1, (d) Level-3 theory for the 1 'S g b ^ 
state at d ~ 7 nm. On the left hand side, the first electron is 
fixed at the center of the bottom dot, while on the right hand 
side, the first electon is fixed at the center of the top dot, as 
indicated by the arrows. 
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FIG. 9: (Color online) Comparison of the DOE calculated 
from (a) Level-1; (b) Level-3 and (c) Level-4 theories for two- 
electron states. In panel (c), both the DOE of the Hubbard 
model (solid lines) and of the Heisenberg model for 1 'Sg state 
(dashed line) are shown. 



E. Comparison of pair correlation functions for 
Levels-1 to 4 theories 



In the previous sections, we compared the energy lev- 
els of two-electron states in several levels of approxima- 
tions to all-bound-state CI results (Level-1). We now 
provide further comparison of Levels 1-4 theories by an- 
alyzing the pair correlation functions and calculating the 
electron-electron entanglement at different levels of ap- 
proximations. 

In Fig. [S] we show the pair correlation functions of 
Eq.@) for the x E^ a) and ^f 1 states at d ~ 7 nm for 
Level-1 and Level-3 theories. The correlation functions 
give the probability of finding the second electron when 



the first electron is fixed at the position shown by the 
arrows at the center of the bottom dot (left hand side of 
Fig. EJ or the top dot (right hand side of Fig. EJ- Level-1 
and Level-2 theories give correlation-induced electron lo- 
calization at large d: for the l Yi^ state, the two electrons 

are localized on different dots, while for the state, 
both electrons are localized on the same doti^* In con- 
trast, Level-3 theory shows delocalized states because of 
the lack of configuration mixing. This problem is shared 
by RHF and LSD approximations. 
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F. Comparison of the degree of entanglement for 
Levels 1-4 theories 

The DOE of the four "£" states are plotted in PigEJfor 
Level-1, Level-3 and Level-4 theories; the DOE of Level-2 
theory are virtually identical to those of Level-1 theory, 
and are therefore not plotted. We see that the Hubbard 
model has generally reasonable agreement with Level-1 
theory while the DOE calculated from Level-3 and Level- 
4 (Heisenberg model) theories deviate significantly from 
the Level-1 theory, which is addressed below: 

(i) The 1 s|, a) state: The Level-1 theory (FigUi), shows 

that the DOE of 1 Sg <1 ' increases with d and approaches 1 
at large d. The Hubbard model of Level-4 theory (Fig|3:) 
gives qualitatively correct DOE for this state except for 
some details. However, Level-3 theory (Fig^) gives 

DOE 5 = because the wavefunction of 1 Yf]f l is a single 
Slater determinant |e) [see Eq.Q]. For the same reason, 
the DOEs of the state in RHF and LSD approxima- 
tions are also zero as a consequence of lack of correlation. 
In contrast, the Heisenberg model of Level-4 theory gives 
5( 1 Sg a ' ) ) = 1. This is because the Heisenberg model as- 
sumes that the both electrons are localized on different 
dots with zero double occupancy, and thus overestimates 
the DOEai^ 

(ii) The state: The Hubbard model gives the 
DOE of the x T,f^ state identical to that of 1 Y, i g a) state. 
This is different from the result of Level-1 theory, espe- 
cially at large inter-dot separations. The difference comes 
from the assumption in the Hubbard model that the en- 
ergy levels and wavefunctions on the top dot and on the 
bottom dot are identical while as discussed in Ref l24l the 
wavefunctions are actually asymmetric due to inhomoge- 
neous strain in the real system. At d > 8 nm, the 
state is the supposition of \E) and \F) configurations in 
the Hubbard model leading to S = 1, while in Level-1 
theory, the two electrons are both localized on the top 
dots (\F)) at d > 9 nm^ resulting in near zero entangle- 
ment. For the same reason discussed in (i), the Level-3 
theory gives l S( 1 Sg 6 ') = 0. 

(hi) The £„ state: Both the Level-3 theory and Hub- 
bard model give iS( 1 E u ) = 1. However, the 5( 1 E M ) of 
the Level-1 theory has more features as the consequence 
of the asymmetry of the system. In contrast to the l T$p 
state, in the 1 S U state, both electrons are localized on 
the bottom dot leading to near zero entanglement at d > 
9 nm. 

(iv) The 3 £ u state: All levels of theories give very close 
results of DOE for qE u state. Actually, in Level-1 theory, 
the DOE of qS u state is only slightly larger than 1, in- 
dicating weak entanglement of the a and tt orbitals (the 
maximum entanglement one can get from the total of six 
orbitals is <S max = log 2 6), while in all other theories (in- 
cluding the Level-2 theory) they are exactly 1 since these 
theories include only two a orbitals. The small coupling 
between a and 7r orbitals is desirable for quantum com- 



putation, which requires the qbits states to be decoupled 
from other states. 



IV. SUMMARY 

We have shown the energy spectrum, pair-correlation 
functions and degree of entanglement of two-electron 
states in self-assembled InAs/GaAs quantum dot 
molecules via all-bound-state configuration interaction 
calculations and compared these quantities in different 
levels of approximations. We find that the correlation be- 
tween electrons in the top and bottom dot is crucial to get 
the qualitative correct results for both the singlet-triplet 
splitting and the degree of entanglement. The single- 
configuration approximation and similar theories such as 
RHF, LSD all suffer from lack of correlation and thus 
give incorrect ground state, singlet-triplet splitting Js~t 
and degree of entanglement. Highly simplified models, 
such as the Hubbard model gives qualitatively correct re- 
sults for the ground state and Js-t, while the Heisenberg 
model only give similar results at large d. These two mod- 
els are written in the dot-centered basis, where the cor- 
relation between the top and bottom dots are carried by 
the single-particle tunneling. However, as a consequence 
of ignoring the asymmetry present in the real system, 
the degree of entanglement calculated from the Hubbard 
model deviates significantly from realistic atomic calcu- 
lations. Moreover the Heisenberg model greatly overes- 
timates the degree of entanglement of the ground state 
as a consequence of further ignoring the electron double 
occupancy in the dot molecule. 
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APPENDIX A: DEGREE OF ENTANGLEMENT 
FOR TWO ELECTRONS 

The entanglement is characterized by the fact that the 
many-particle wavefunctions can not be factorized as a 
direct product of single-particle wavefunctions. An en- 
tangled system displays non-locality which is one of the 
properties that distinguishes it from classic systems. So 
far, the only well established theory of entanglement per- 
tains to two distinguishable particles, 32 ' 34 (e.g. electron 
and hole). For a system made of two distinguishable par- 
ticles (A,B), the entanglement can be quantified by von 
Neumann entropy of the partial density matrix of cither 
A or jj, 31 i 32 i 33 

S(A, B) = -Tr (p A log 2 p A ) = -Tr (p B log 2 p B ) , (Al) 
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where S(A, B) is the DOE of the state, pa and ps are the 
reduced density matrices for subsystems A and B. An 
alternative way to define the DOE for two distinguishable 
particles is through a Schmidt decomposition, where two- 
non-identical-particle wavefunctions are written in an bi- 
orthogonal basis, 

*(A,B)=J2*i\iA)®\iB), (A2) 

i 

with Ai > and ^ A? = 1 . The number of nonzero Xi 
is called the Schmidt rank. For a pure state *i>(A,B) of 
the composite system (A, B), we have, 

PA = J2 X * i^x^i 

i 

PB = ^A 2 \iB)(i B \. (A3) 

i 

It is easy to show from Ea. lAll that the DOE for the two 
distinguishable particles is, 

S(A,5) = -^A 2 log 2 A 2 . (A4) 

i 

We see from Ea. (|A2|) that when and only when the 
Schmidt rank equals 1, the two-particle wavefunction can 
be written as a direct product of two single-particle wave- 
functions. In this case, we have A = 1, and S(A, B) = 
from Ea. (|X4l) . 

A direct generalization of DOE of Ea. (|A4|) for two iden- 
tical particles is problematic. Indeed, there is no general 
way to define the subsystem A and B for two identi- 
cal particles. More seriously, since two-particle wave- 
functions for identical particles are non-factorable due 
to their built-in symmetry, one may tend to believe that 
all two identical fermions (or Bosons) are in entangled 
Bell stated 2 However, inconsistency comes up in the lim- 
iting cases. For example, suppose that two electrons are 
localized on each of the two sites A and B that are far 
apart, where the two electrons can be treated as distin- 
guishable particles by assigning A and B to each electron, 
respectively. A pure state ^ that has the spin up for A 
electron and spin down for B electron is ^(x^x^ — 
1/v / 2[0 j 4t(xi)</>b|(x 2 ) - 0at(x2)</>bx(xi)]. At first sight, 
because of the anti-symmetrization, it would seem that 
the two electron states can not be written as a direct 
product of two single particle wavefunctions, so this state 
is maximally entangled. However, when the overlap be- 
tween two wavefunctions is negligible, we can treat these 
two particles as if they were distinguishable particles and 
ignore the anti-symmetrization without any physical ef- 
fect, i.e. \&(xi,X2) = j 4t( x i)</ ) b|(x2). In this case, 
apparently the two electrons are unentangled. More in- 
triguingly, in quantum theory, all fermions have to be 
anti-symmetrized even for non-identical particles, which 
does not mean they are entangled. 

To solve this obvious inconsistency, alternative mea- 
sures of the DOE of two fermions have been proposed and 



discussed recently) 35 i 36 i 37 i 38 i 39 i °i but no general solu- 
tion has been widely accepted as yet. Schliemann et a.' 35 
proposed using Slater decomposition to characterize the 
entanglement (or, the so called "quantum correlation" in 
Ref. 13511 of two fermions as a counterpart of the Schmidt 
decomposition for distinguishable particles. Generally a 
two-particle wavefunction can be written as, 

* = 5^w ii |i)®|i), (A5) 

i,3 

where \i), \j) are the single particle orbitals. The coeffi- 
cient u>ij must be antisymmetric for two fermions. It has 
been shown in Ref. I35ll42l that one can do a Slater de- 
composition of u>ij similar to the Schmidt decomposition 
for two non-identical particles. It has been shown that uj 
can be block diagonalized through a unitary rotation of 
the single particle statesj2£4£ i.e., 

J = Uuj C/t = diag[Zi, Z 2 , ■ • • , Z r , Z ] , (A6) 

where, 

Zi = ( X o ) ■ < A7 > 

and Zq = 0. Furthermore, ^i^f—l, and z% is a non- 
negative real number. A more concise way to write down 
the state ^ is to use the second quantization representa- 
tion, 

* = I>/Ll/£|0>, (A8) 

i 

where, / 2 j_i and /i are the creation operators for modes 
2i — 1 and 2i. Following Ref|42j, it is easy to prove that 
zf are eigenvalues of upuj. The number of non-zero Zi is 
called Slater ranki^ It has been argued in Refl35l that if 
the wavefunction can be written as single Slater determi- 
nant, i.e., the Slater rank equals 1, the so called quantum 
correlation of the state is zero. The quantum correlation 
function defined in Ref |3^ has similar properties, but nev- 
ertheless is inequivalent to the usual definition of DOE. 

Here, we propose a generalization of the DOE of 
Eci. (|A4(l to two fermions, using the Slater decomposi- 
tions, 

S = -5> 2 log 2 z 2 . (A9) 

i 

The DOE measure of Eq. I|A9|) has the following proper- 
ties: 

(i) This DOE measure is similar to the one proposed 
by Paskauskas et aim* and Li et aly 9 except that a differ- 
ent normalization condition is used. In our approach, the 
state of Slater rank 1 is unentangled, i.e., 5=0. In con- 
trast, Paskauskas et al^ and Li et al^> concluded that 
the unentangled state has S = In 2, which is contradic- 
tory to the fact that for distinguishable particles, an un- 
entangled state must has 5=0. In our approach, the max- 
imum entanglement that a state can have is S = log 2 N, 
where N is the number of single particle states. 
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(ii)The DOE measure of Eq. i|A9|l is invariant under 
any unitary transition of the single particle orbitals. Sup- 
pose there is coefficient matrix lo, a unitary transforma- 
tion of the single particle basis leads to a new matrix 
u>' = UojW and u'*oj' = UutujW. Obviously, this 
transformation would not change the eigenvalue of uj^uj, 
i.e., would not change the entanglement of the system. 

(hi) The DOE of Eq. i|A9(l for two fermions reduces to 
usual DOE measure of Eq. (|A4|) for two distinguishable 
particles in the cases of zero double occupation of same 
site (while the DOE measure proposed by Paskauskas 
et alr^ and Li et al.— does not). This can be shown 
as follows: since the DOE of measure Eq. (|A9|) is ba- 
sis independent, we can choose dot-localized basis set, 
(which in the case here is the top (T) and bottom (B) 
dots [Fig^b)]), such that the antisymmetric to matrix in 
the dot-localized basis has four blocks, 



where, ipi is the i-th molecular orbital, and IA is a unitary 
matrix, i.e. WIA = I. To obtain a set of well localized 
orbitals, we require that the unitary matrix IA maximizes 
the total orbital self-Coulomb energy. The orbitals ful- 
filling the requirement are approximately invariant un- 
der the changes due to coupling between the dots^i For 
a given unitary matrix IA, the Coulomb integrals in the 
rotated basis are, 



f"7i>»72 



E 



u 



■11 



Vi ,1 1 k,l ' 



(B2) 



where TVi are the Coulomb integrals in the molecular ba- 
sis. Thus, the total self-Coulomb energy for the orbitals 
{X^} is: 



_ j UlTT —^TB 
OJTB LtlBB 



(AlO) 



T) i,j,k,l 



>1J 



(B3) 



where, lott is the coefficient matrix of two electrons both 
occupying the top dot, etc. If the double occupation is 
zero, i.e., two electrons are always on different dots, we 
have matrices ujtt — ^bb—0. It is easy to show that 
uj* lo has two identical sets of eigenvalues zf, each are 
the eigenvalues of u)^ bt usbt- On the other hand, if we 
treat the two electrons as distinguishable particles, and 
ignore the anti-symmetrization in the two-particle wave- 



functions, we have ps 



,t 

TB 



u>tb and px 



BT 



lubt- 



It is straightforward to show that in this case Eq. 
and Eq. (|A4|) are equivalent. 



APPENDIX B: CONSTRUCTION OF 
DOT-CENTERED ORBITALS 



The procedure of finding the unitary matrix IA that max- 
imizes t/tot is similar to the procedure given in Ref. 
where the maximally localized Wannier functions for ex- 
tended systems are constructed using a different criteria. 
Starting from IA = I, we find a new IA = I + Se that 
increases E/tot- To keep the new matrix unitary, we re- 
quire 5e to be a small anti-Hermitian matrix. It is easy 
to prove that 



SU U 



Se 



3,1 



(B4) 



When we solve the single-particle Eq.Q for the QDM, 
we get a set of molecular orbitals. However sometimes 
we need to discuss the physics in a dot-localized basis 
set. The dot-localized orbitals Xn can be obtained from 
a unitary rotation of molecular orbitals, 



N 



Xn 



(Bl) 



G* i . By choosing Sei 



and to verify that Gjj =- . . ( , 
—eGij, where e is a small real number, we always have 
(to the first-order of approximation) AZ7 t ot = e||G|| > 0, 
i.e. the procedure always increases the total self-Coulomb 
energy. To keep the strict unitary character of the IA 
matrices in the procedure, the IA matrices are actually 
updated as IA «— IA exp(— eG), until the localization is 
achieved. 
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